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The angle-dependences of the magnetoresistance of two different isotopic substitutions (deuterated 
and undeuterated) of the layered organic superconductor k-(ET)2Cu(NCS)2 are presented. The 
angle dependent magnetoresistance oscillations (AMRO) arising from the quasi-one-dimensional 
(Q1D) and quasi-two-dimensional (Q2D) Fermi surfaces in this material are often confused. By 
using the Boltzman transport equation extensive simulations of the AMRO are made that reveal the 
subtle differences between the different species of oscillation. No significant differences are observed 
in the electronic parameters derived from quantum oscillations and AMRO for the two isotopic 
substitutions. The interlayer transfer integrals are determined for both isotopic substitutions and a 
slight difference is observed which may account for the negative isotope effect previously reported pj. 
The success of the semi-classical simulations suggests that non-Fermi liquid effects are not required 
to explain the interlayer-transport in this system. 
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I. INTRODUCTION 



k-(ET)2Cu(NCS)2 is probably the most popular and 
best characterised material out of all the organic charge 
transfer salts based on the ET molecule. Its attrac- 
tion to experimentalists lies in its exceedingly simple 
Fermi surface, which consists of two elliptical quasi-two- 
dimensional (Q2D) pockets and apair of warped quasi- 
one-dimensional (Q1D) sheets 0,0 (see Figure 0. The 
prospect of understanding the complex transport proper- 
ties of the organic salts seems more within reach for this 
material than for others that show similar behaviour but 
have more complicated Fermi surfaces. 

Several theoretical models of the unconventional su- 
perconductivity observed in /c-(ET)2Cu(NCS)2 and re- 
lated materials suggest that the superconducting pairing 
mechanism may be mediated by antiferromagnetic spin 
fluctuations 0, M E| . These models are found to be sen- 
sitive to the degree to which the Fermi surface of the 
material can nest; the higher the nestability the more 
likely this pairing is to be successful. Two-dimensional 
Fermi surfaces are clearly better able to nest than three- 
dimensional ones, and so tests of the dimensionality of 
k-(ET)2Cu(NCS)2 also test these theoretical models. 

In this paper the low temperature angle-dependence 
of the magnetoresistance in deuterated and undeuter- 
ated samples of k-(ET) 2 Cu(NCS)2 is studied in detail in 
magnetic fields significantly higher than in-plane upper 
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critical field. This is the first time that comprehensive 
measurements like these have been made. Their purpose 
is to completely determine the parameters that define the 
transport in this material, to locate any differences be- 
tween these parameters for the two isotopic substitutions 
that might shed light on the disparity between their su- 
perconducting critical temperatures an d to address 
the question of whether it is possible to describe all as- 
pects of the normal state transport within the bounds of 
Fermi liquid theory. 

As with all the organic conductors in this class, the 
ET molecules form the highly conducting layers, sepa- 
rated by layers of the anion, with the long axis of the ET 
molecule at a small angle to the interlayer direction. In 
the K-phase salts the ET molecules associate into pairs, 
or dimers, each of which collectively donates one electron 
to the anions, leaving behind a mobile hole pj. There are 
two dimers, and thus two holes per unit cell, and so, be- 
cause the dispersion is nearly isotropic in the fee-plane, 
this leads to a roughly circular Fermi surface which has 
the same area as the first Brillouin zone Q. The Bril- 
louin zone itself reflects the rectangular cross-section of 
the unit cell and the Fermi surface cuts the Brillouin 
zone boundaries on its long side. At these points a gap 
opens up which splits the Fermi surface into the Q1D and 
Q2D section 
of Figure ^ 

The shape of the Fermi surface in the fc^fcy-plane has 
been confirmed by the observation of magnetic quantum 
oscillations 0, H|. The frequency of the quasiparticle 
orbits about the circumference of the Q2D pockets (a- 
orbits) is found to be 600 T which corresponds to about 
15 % of Abz, the area of the cross-section of the first 
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FIG. 1: Top: Cross-section of the Fermi surface of k- 
(ET)2Cu(NCS)2 in the k x k y -p\a.ne. The shape is defined by 
three transfer integrals in the highly conducting layers (see 
Equation 0. Bottom: The same Fermi surface in three di- 
mensions. The interlayer warping is defined by the transfer 
integral t a (see Equation [(J, and is exaggerated for clarity. 

Brillouin zone in the k x k y plane Above 20 T mag- 
netic breakdown is observed as it becomes possible for 
some of the quasiparticles to bridge the energy gap be- 
tween the two Fermi surfaces and make the large /3-orbit, 
whose frequency is found to be 3920 T This corre- 
sponds to an area equal to Abz to within a few percent, 
as would be expected from the considerations discussed 
above. 

k-(ET)2Cu(NCS)2 has a monoclinic structure and the 
transfer integral between the layers lies parallel to the a 
lattice parameter which is inclined at angle of 110.3° to 
the highly conducting frc-planes 0. This transfer inte- 
gral, i a , is much smaller than those within the planes and 
results in a slight warping of the Fermi surface perpen- 
dicular to the direction of t a in /c-space. This is shown 



in the lower portion of Figure ^ The validity of this pic- 
ture of the Fermi surface as a three-dimensional object is 
discussed in Reference 0. 



II. OVERVIEW OF THE FEATURES 
OBSERVED IN MAGNETOTRANSPORT 

Figure [3 shows two typical ^-dependences (where 8 is 
the angle between the magnetic field and the normal 
to the conducting layers) of the interlayer magnetore- 
sistance of k-(ET) 2 Cu(NCS)2 in fixed magnetic fields of 
27 T and 42 T and at an azimuthal angle of 149°. In such 
high fields a whole host of features are observed in an 
interlayer transport measurement of k-(ET)2Cu(NCS)2, 
e.g. Shubnikov de Haas (SdH) oscillations, magnetic 
breakdown, and Q1D and Q2D angle dependent mag- 
netoresistance oscillations (AMRO). This means that for 
a typical 0-rotation the magnetoresistance is rich in fea- 
tures as Figure |3 illustrates. 

The upper plot shows the data taken at 42 T and the 
lower at 27 T. In the upper plot the field perpendicular 
to the layers around 6 = 0° is sufficient for the effects of 
magnetic breakdown to be observed and the fast SdH os- 
cillations due to the /3-orbit are clearly seen. The slower 
oscillations due to the a-orbit are seen in both plots and 
persist to higher angles. The amplitudes of these os- 
cillations are modulated and they disappear at certain 
6>-angles, these nodes are known as spin-zeroes and are 
caused by Zeeman splitting of the Landau levels ^(| . An 
analysis of this effect is dealt with in Section IV Bl 

The positions of the features at 0-angles greater than 
about ±70° are seen to be independent of the magni- 
tude of the magnetic field, which reveals them to be 
AMRO of one variety or another. Four different types of 
AMRO are possible in the interlayer resistivity (p zz ) of 
k-(ET) 2 Cu(NCS) 2 . These are the Danner-Kang-Chaikin 
oscillations [ll| , the third angle effect 0] an d the Lebed 
magic angle effect [l3j which all arise from orbits on 
the Q1D Fermi surface section, and the Yamaji oscil- 
lations arising from orbits on the Q2D Fermi surface sec- 
tion 

mull- in th ° 

semi-classical picture all the AMRO 
are caused by the degree to which the velocity compo- 
nents of the quasiparticle are averaged over the series of 
orbits that appear at a certain inclination angle. In par- 
ticular, the orbits that are possible in the region of the 
Yamaji angles are very successful in averaging the inter- 
layer velocity towards zero, thus peaks are seen in the 
interlayer resistance [l4lll5j. In contrast, the orbits that 
occur at the Lebed magic angles are not as successful at 
averaging the interlayer velocity towards zero as those 
possible at the other angles and so dips in p zz are ob- 
served 0] . There are other theories that can explain the 
effects observed at the Lebed magic angles. Lebed's own 
argument describes electron-electron correlations whose 
magnitudes change when the field is directed along the 
magic angles 0] . Another theory has regions of /c-space 
where the scattering rate takes a large value (Fermi sur- 
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FIG. 2: Typical ^-dependence of the magneto-resistance of «-(ET)2Cu(NCS)2. The data shown is for a hydrogenated sample at 
490 mK, 4> — 149°(where <j> is the azimuthal angle), 27 T (lower) and 42 T (upper). The data have been offset for clarity. Some 
representative features are indicated; SdH oscillations due to the Q2D pockets (a) and the breakdown orbit (/3); spin-zeroes 
in the SdH amplitudes (SZ); the onset of the superconducting transition (SC); angle-dependent magnetoresistance oscillations 
(AMRO), whose positions are field independent; and the resistive peak in the presence of an exactly in-plane magnetic field 
(In- plane Peak). The inset diagram is included to illustrate the measurement geometry. 



face hotspots) accounting for the AMRO ^t|- However, 
such theories are complicated and need only be invoked 
when the semiclassical approach fails to account for the 
experimental observations. It will be shown by the sim- 
ulations described in Section IV CI that the semi-classical 
explanation is sufficient in the case of k-(ET)2Cu(NCS)2- 

In the upper plot of Figure [21 a small peak is observed 
when the field lies very close to the in-plane direction, 
9 ss 90°. This is the in-plane peak feature mentioned in 
Reference @. It will be discussed further in Section IVD1 
Around 9 — 90° in the lower plot the in-plane peak is 
obscured by the large dip that indicates the onset of a 
superconducting transition. This occurs because there is 
a considerable anisotropy in the upper critical field of this 
material, and a field of 27 T is not sufficient to suppress 
the superconducting state when applied in a nearly in- 
plane direction |l8j. 

It should also be noted from Figure [3 that the ampli- 
tude of all the features in the magnetoresistance increase 
with increasing field, and that the plots are not symmet- 



rical about 9 = 0°, reflecting the monoclinic symmetry 
of the crystal structure. 



III. PARAMETERISING THE FERMI SURFACE 

It has been shown that the measured intralayer Fermi 
surface of k-(ET)2Cu(NCS)2 can be reproduced using a 
dispersion relation derived from a tight binding model 
using the ET dimer as its base unit 0, S 13- I n this 
way the intradimer transfer integral, td, can be ignored 
and shape of the Fermi surface depends upon interdimer 
transfer integrals i&, t c i, i C 2 and the Fermi energy, E-p. 
The dispersion found in this manner is known as the ef- 
fective dimer model and is given by 



E(k) = 2t b cos(k x b) 
k x b 
~2 



± cos I — Jt 2 cl + t 2 c2 + 2t cl t c2 cos(k y c) (1) 
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where the + and — signs result in the Q1D and the Q2D 
sections of the Fermi surface respectively 0, § . 

The effective dimer model is used in the semi-classical 
calculations of Sections IV CI and IV Dl These do not take 
account of quantum effects such as Shubnikov de Haas 
oscillations or magnetic breakdown, and are assumed to 
be in the "low-field" region where breakdown does not 
occur. It will be seen that this is a reasonable assump- 
tion in both cases. This means that the effect of the 
energy gap, i.e. the difference between t c \ and t C 2, can 
at first be neglected and the shape of the in-plane Fermi 
surface depends only upon the ratios Ev/t c and %/t c , 
where t c is an average of t c \ and t c2 0, 13 • 

It is possible to obtain values for these ratios by ad- 
justing them to reproduce the areas of the a and (3 Fermi 
surface orbits. Once this is done tt, t c and E-p can be 
uniquely specified by fitting to the effective mass of the 
/3-orbit as found from SdH oscillations, using the expres- 
sion SAp/SE — 2Trm*p/H 2 Note that it also possible 
to fit to the mass of the a-orbit to obtain slightly dif- 
ferent results. However, as the masses are derived from 
quantum oscillations, they are orbitally averaged, and so 
the a-mass will be dominated by the extremely pointed 
regions of the Q2D pockets. The breakdown orbit does 
not have these pointed regions and thus it is the /3-mass 
that is used in the fitting procedure. 

The energy gap can now be reintroduced in order to 
specify t c \ and £ c2 . In the region of the gap, cos(fc c c) = 
— 1 and cos(fcb&/2) w 0.5, so that E g w 2(t c i — i c2 ). 
From the magnetic breakdown, E g is estimated to be 
7.8 meV [9l |20|. Any inaccuracies in this value will not 
lead to errors in the size of the Fermi surface produced 
by Equation ^ but could lead to small discrepancies in 
the exact dimensions of the a-pocket. 

The Fermi surface parameters obtained in this man- 
ner are as follows: tf, = 14.87 meV; t c \ — 26.65 meV; 
t c2 = 22.75 meV; and E F = -19.12 meV. Note that 
these values of t are effective transfer integrals, and in- 
corporate the effects of electron-phonon and electron- 
electron interactions as they are derived from magnetic 
quantum oscillation data (see [5l| for a discussion) . Note 
also that E-p is taken relative to the zero energy of the 
effective dimer model and not the bottom of the band. 
It should therefore not be quoted as the Fermi energy of 
k-(ET) 2 Cu(NCS) 2 . 



IV. EXPERIMENTAL DETAILS 

Four single crystal samples of k-(ET) 2 Cu(NCS) 2 were 
used in this study, made using an electrocrystallisation 
method 0, and references therein] . All of the samples are 
black platelets of the order of 0.7x0.5x0.1 mm 3 , with the 
plane of the plate corresponding to the highly conducting 
layers. In one of the samples the eight terminal hydrogens 
of the ET molecules were substituted by deuterium. In 
what follows the deuterated crystal will be referred to as 
d8, and the hydrogenated crystals as h8. 



The magnetoresistance measurements were made using 
standard 4-wire A.C. techniques (/ = 50 — 180 Hz) with 
the current applied in the interplane direction (7=1 — 
20 (j,A). 

All the samples were mounted on a two-axis rotator in 
a 3 He cryostat. In this rotator it is possible to continu- 
ously change the 0-angle, the angle between the magnetic 
field and the highly conducting frc-planes, and discretely 
change the plane of rotation, described by the azimuthal 
angle, <p. An angular calibration technique similar to that 
described in Reference was used when misalignments 
of the sample that occur during cooling were found to be 
significant. Temperatures down to 0.5 K are readily ac- 
cessible. 



V. RESULTS AND DISCUSSION 
A. The Shubnikov de Haas oscillations 

Using the results of a fast Fourier transform analysis of 
several SdH measurements the fundamental frequencies 
were found to be F a {h8) = 599 ± 3 T; F p (h8) = 3860 ± 
6 T; F a {d8) = 598 ± 3 T; and 7>(d8) = 3871 ± 10 T, 
all of which are in reasonable agreement with previous 
results dlH. 

The a-mass of h8 k-(ET) 2 Cu(NCS) 2 has previously 
been found to be m* a (h8) = 3.5 ± 0.1 m e Q using 
a Lifshitz-Kosevich analysis of the temperature depen- 
dence of the SdH amplitudes Using the same 
method the equivalent c?8-mass is to* (g?8) = 3.6 ±0.1 me. 

The Lifshitz-Kosevich analysis can also be applied to 
the field-dependence of the SdH amplitudes at a con- 
stant temperature to find values for the scattering time 
and the breakdown field. Using the technique outlined in 
Reference [20] (but correcting the erroneous minus sign 
that prefixes the breakdown term in that reference) the 
amplitudes, A, are fitted with the function 
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: sinh( 7j r/7j) 
IjT/B 



D 
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(2^ix*hk B )/(he), p 2 



(2) 

ex-p(—B j/B) and 



where 7^ 

q 2 = 1 — p 2 [lj| . Tb is the Dingle temperature and is 
proportional to the scattering rate, fi* = m*/m c , m c 
is the mass of an electron, lj is the harmonic index of 
the orbit, n\j is the number of magnetic breakthrough 
points on the orbit, and n 2j - is the number of Bragg reflec- 
tion points 0|- Note that this is derived from the two- 
dimensional form of the Lifshitz-Kosevich formula with 
the magnetic field directed perpendicular to the highly 
conducting layers. The amplitudes are obtained from the 
fast Fourier transform spectrum of the oscillating part of 
the resistance. The field window over which the Fourier 
transform is performed specifies the value of B in the 
above equation such that B~ x = (75-j -1 + B 2 ~ 1 )/2, where 
Bi and 7? 2 are respectively the start and end points of 
the field window 1241. 
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The functional form of the Tb and Bq terms in the 
Lifshitz-Kosevich formula are similar and so to obtain 
a satisfactory fit the amplitudes of the first and second 
harmonics of the a-frequency must be fitted simultane- 
ously 20]. In this case n\{a) = n\{2a) = 0, 712(a) = 2 
and U2(2a) — 4. Thus, at high fields these amplitudes 
are attenuated as quasiparticles are able to tunnel across 
the gap to the /3-orbit. The fits are shown in Figure |3| 
and the values obtained are r(h8) — 2.3 ±0.2 ps; r(d8) = 
2.4 ±0.2 ps; B (h8) = 58 ±9 T; and B (d8) = 39 ±10 T. 
The values of the scattering time obtained from the high 
field fits are in close agreement with those obtained from 
a fit to the low-field data where the effects of breakdown 
may be neglected. 

Note that all the results for the d8 sample are the same 
as those for the h8 to within the experimental errors. 
The large errors on the values of the breakdown field 
and the discrepancies between these values and those 
independently obtained from similar data (Bo(h8) = 
41 ±7 T [H ) and even the same data (B (h8) = 41 ±5 T 
and B (d8) = 30 ±5 T jH) serve to highlight the limita- 
tions of the Lifshitz-Kosevich formula at high fields [25| , 
where not only is there competition between two func- 
tionally similar terms, but also the amplitudes of the 
quantum oscillations become very large. 



B. The effect of spin-splitting 

The energy levels of a quasiparticle in a metallic sys- 
tem subjected to an applied magnetic field are defined 
by Landau quantisation and the Zeeman effect, and are 
given by 



E 



, 1 . heB cos 9 , 1 4 _, 
( n +n) ; ±7J>bB, 



(3) 



where n is the Landau level index and g* is the effective 
g-factor [2(j. Increasing the 0-angle reduces the separa- 
tion between Landau levels by reducing the field perpen- 
dicular to the highly conducting planes, Bcos9. When 
B cos 9 is such that the spin- up and spin-down sections 
of different Landau levels are degenerate then the sepa- 
ration between successive energy levels is equal to huo c . 
At this angle the SdH oscillations having the fundamen- 
tal frequency, F, will dominate, taking their maximum 
amplitude. However, when B cos 9 is such that the spin- 
up and spin-down sections of different Landau levels are 
equally spaced at ^fko c , then the dominant oscillations 
will be those with frequency 2F and the amplitude of 
the fundamental oscillations will be at a minimum |10| . 
These two situations arc known as spin-maxima and spin- 
zeroes respectively. 

It is easy to show that the conditions for spin-zeroes 
and spin maxima are given by |10| 



g*HBB = j- 
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FIG. 3: The result of fitting the Fourier amplitude, A, of 
the h8 (top) and d8 (bottom) a-frequency (squares) and its 
second harmonic (circles) to the two-dimensional Lifshitz- 
Kosevich formula at constant temperature, over a field range 
of 12 - 44 T, using the technique outlined in the text and 
Reference [2(j. This kind of analysis yields values for the 
scattering time and the magnetic breakdown field. The in- 
sets show the result of fitting the temperature dependence 
of the h8 (top) and d8 (bottom) a-frequency amplitude over 
a constant (low-)field interval. The h8 and d8 data at each 
temperature were taken simultaneously and the fact that the 
insets are similar is an indication that the effective masses of 
the two isotopic substitutions are also similar. 



As has already been mentioned, this effect can be ob- 
served as a modulation of the a-frequency SdH oscilla- 
tions when a crystal of k-(ET)2Cu(NCS)2 is rotated in 
a fixed field. Figure 01 shows a typical section of such a 
rotation for an ft.8-sample, and several spin-zero angles, 
9j, are marked with arrows. The inset shows a plot of 
(cos^j) -1 versus j-index. This dataset is a summary of 
a large number of spin-zero positions measured at many 
different values of the azimuthal angle, at fields of 27 and 
42 T, and in two different single crystals. Using the gra- 
dient of the straight line fit shown, the product of the 
effective g-factor and //*, where \x* — m* /m e , is found to 
be g*nl(h8) = 5.22 ± 0.56. 
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FIG. 4: The spin-zero effect as seen in a typical ^-dependence 
of the magnetoresistance of k-(ET)2Cu(NCS)2- The data 
shown is for an h8 sample at 590 mK, (f> — 90.8°, and 27 T. 
The spin-zero angles, 6j are indicated by arrows. The data 
points in the inset are a summary of the 9j observed at a num- 
ber of different azimuthal angles, magnetic fields, and in two 
different single crystal samples. The broken line is a straight 
line fit to these data. 



The dataset for d8 is not as extensive as that for 
h8, nevertheless a good fit is still achieved, yielding 
g*fj,* a (d8) = 5.24 ± 0.65. The values for hB and d8 are 
identical within the errors. 

The value for (h8) obtained here app ears to be in 
good agreement with that of Reference [27j , obtained by 
fitting three spin-zero points from de Haas-van Alphen 
data. In that reference the authors assume that g* = 2 
and that the mass obtained from the spin-zero effect is 
renormalised by electron-electron interactions, but not 
electron-phonon interactions, they then use the differ- 
ence between this mass and that derived from a Lifshitz- 
Kosevich analysis of quantum oscillations to specify the 
electron-phonon coupling constant [l]}. However, as the 
g* /i* values obtained in this manner may not be renor- 
malised in the same way as the effective masses found 
from the thermodynamic variation of the quantum oscil- 
lation amplitudes or the g-factors obtained from electron 
spin resonance it is not advisable to separate g* and fi* 
in this fashion. 

It is possible, using the experimentally determined 
value of to make a comparison of the spin and 

Landau level splittings. The ratio of the splittings at 
9 = 0° is given by, 



g*/d B B _ g*n* 



hljj r 



(5) 



using uo c — eB/m* and /iB = eh/2m e . And so for both 
h8 and d8 the ratio of the spin to Landau level split- 
ting is around 2.6, which results in the energy level di- 
agram shown in Figure It is seen that as well as the 
energy splitting of ftw c between one spin split level and 
its equivalent from the next highest Landau level, there 
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FIG. 5: The energy level spectrum of k-(ET) 2 Cu(NCS) 2 at 
6 = 0° as deduced from g* fj,^ ~ 5.2, measured at ~ 0.5 K 
in the low-field region. There are two splittings, one of h(jj c , 
which results in the SdH oscillations of the fundamental fre- 
quency, and one of 0.4 hw c , which at sufficiently low tempera- 
tures and high fields will result in the observation of harmonics 
of the fundamental frequency. 



exists another splitting, of 0.4 huj c , arising from the dif- 
ference between the spin-up of one Landau level and the 
spin-down of the Landau level three places up. 



The angle-dependent magnetoresistance 
oscillations 



1. Boltzmann transport simulations 

An analysis of the angular effects in k- 
(ET) 2 Cu(NCS)2 is complicated by the co-existence 
of Q1D and Q2D Fermi surfaces. The method by which 
the different types of AMRO either dominate or super- 
pose over one another is not at all clear, depending as it 
does on unknowns such as the relative effective masses 
and carrier densities of the quasiparticles on the Q1D 
and Q2D sections. To further complicate an AMRO 
investigation it should be noted that in general Lebed 
magic angles and Yamaji oscillations can be analysed 
in very similar ways. For example, if the resistive peaks 
that lie between dips caused by the Q1D Lebed magic 
angle effect are accident ly mistaken for Q2D Yamaji 
oscillations it is possible, as will be shown later, to 
obtain the dimensions of a closed Fermi surface pocket 
that may appear reasonable, but is incorrect. For this 
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reason, when measuring samples whose Fermi sur- 
face is uncharted, the Lebed magic angles and Yamaji 
oscillations are best used in conjunction with other 
Fermi surface effects such as Danner-Kang-Chaikin, 
or quantum oscillations, which specify exclusively the 
nature of the Fermi surface from which they arise. 

In the sample under review here, the presence of both 
Q1D and Q2D sections of Fermi surface is not in question 
as it is demonstrated convincingly by the magnetic break- 
down observed in the SdH effect. However, in order to 
make sense of the AMRO data measured experimentally, 
some method of separating the oscillations arising from 
the two sections is required. This is achieved by making 
detailed, semi-classical simulations of the interplane re- 
sistivity resulting from the Q1D and Q2D Fermi surfaces. 
A suite of programs were therefore developed which used 
Fortran (for operational speed) to solve the equations 
of motion for any specified Fermi surface and field ori- 
entation and use the results of this to find a numerical 
solution to the Chambers formula (Equation [7| below). 
This software was applied to model Fermi surfaces and 
the AMRO results were seen to agree with theoretical 
predictions. 

It is necessary here to simulate the angle-dependent 
effects observed in k-(ET) 2 Cu(NCS) 2 . To this end an 
equation that describes the entire Fermi surface of this 
material throughout the first Brillouin zone is formu- 
lated; 

E(k) = 2t b cos(k x b) 



This argument is formalised in the isothermal solu- 
tion to the Boltzmann transport equation known as the 
Chambers formula; 



± cos ( ^ ) \jt 2 cl + t 2 c2 + 2t c it c2 cos(kyc) 



-2t a cos 



7T 7T 

k z acos(/3 - -) - k y asm(/3 - -) (6) 



in which the conducting layers of the (ET) molecules lie 
in the be- (or xy-) plane, f3 is the angle between the crys- 
tallographic a and c directions, and the z-axis lies along 
the interlayer direction. The first two terms of the equa- 
tion are simply the effective dimer model that was dis- 
cussed previously and which is known to accurately de- 
scribe the intralayer Fermi surface. The last term is a 
tight binding representation of the interlayer dispersion. 
For the purposes of the simulations the values of ti,, t c i, 
t C 2 and the Fermi energy are set to be those quoted ear- 
lier, derived from a consideration of experimental results. 
The transfer integral t a is set to be 0.04 meV. This is the 
value resulting from a preliminary investigation of the 
in-plane peak effect 0, |28j . 

As a quasiparticle moves across the Fermi surface un- 
der the influence of the magnetic field its component of 
velocity in a given direction will vary as it negotiates the 
various Fermi surface contours and corrugations (its total 
velocity remaining perpendicular to the Fermi surface at 
all times) . Considering an entire Fermi surface filled with 
orbiting quasiparticles it is not difficult to see that the 
combination of these varying velocity components leads 
to a conductivity that depends strongly on the nature of 
the orbits and the geometry of the Fermi surface. 



e 

47T 3 



dk 3 {- 
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de 



)«i(k,0) / ty(k, t)e*/ r di, (7) 



where cry is a component of the conductivity tensor, /o is 
the unperturbed quasiparticle (Fermi-Dirac) distribution 
function, vt and Vj are velocity components and 1 /r is the 
k-independent scattering rate |29j |. This equation repre- 
sents a velocity-velocity correlation function between the 
i th component of the initial velocity, «i(k, 0), integrated 
over all possible starting points on the Fermi surface, and 
Vj(k, t), the j th component of the velocity of a quasipar- 
ticle averaged over the duration of its orbit. The expo- 
nential term represents the probability of a quasiparticle 
scattering from its trajectory so that it no longer con- 
tributes to the conductivity. 

Armed with the Chambers formula and Equation 
it is now possible to relate the way in which the pro- 
gram that simulates the interplane resistivity proceeds 
(a discussion of the possible errors that might creep in 
is left until the end): first a quasiparticle is placed at 
point on the Fermi surface and its velocity components 
are found and recorded by differentiating Equation^] ac- 
cording to v = ftr 1 Vk-E(k) |3*flj . Next, the Lorentz force 
(hdk/dt = — ev x B [3(J) for a given inclination of the 
magnetic field is allowed to act upon the quasiparticle for 
a short time so that it moves to a new position on the 
Fermi surface. Here its velocity components are again 
recorded and the process is repeated a large number of 
times so that a Fermi surface orbit is mapped out. 

The time- integral in the Chambers formula is obtained 
by multiplying each value of each velocity component by 
exp(— t/r) and adding the like components together. The 
scattering time, r, is chosen to be 3 ps in order to re- 
flect that measured from quantum oscillations and high 
frequency conductivity measurements .9] and the whole 
orbit is recorded over a time t — 8r, by which point 
more than 99.96% of the quasiparticles have been scat- 
tered. The time interval, 5t, between points on an orbit 
was set to be quite small (never larger than 0.002t) so 
that the sum of velocity components might approximate 
well to the integral in the Chambers formula. The pro- 
gram returns to the correct position any orbit that has 
a tendency to move off the Fermi surface using a sub- 
routine that exploits the Runge-Kutta method of solving 
ordinary differential equations. 

The time-averaged velocity components are now mul- 
tiplied by the relevant velocity component from the start 
of the process, i.e. t — 0, and weighted by the density 
of states and the Fermi surface area represented by the 
orbit. This routine is repeated for a large grid of starting 
points that span the entire first Brillouin zone, and the 
results are summed. In this way the integral over the 
Fermi surface in the Chambers formula is accomplished, 
and each component of the conductivity tensor is calcu- 
lated. The results are combined to yield the interplane 
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resistivity by inverting the conductivity tensor; 



"xx^yy ^xy^yx ^ 

&xx&yy&zz ^xx^yz® zy "h @xz&yx@zy &xy&yz&zz ^xy^yz^xz @ xz@yy@ zx 



This method can be used to calculate the resistivity 
at any value of 9, 4> and B. The Fermi surface resolu- 
tion chosen, i.e. the number of orbits sampled, must be a 
compromise between the accuracy of the results and the 
speed of calculation. It is found that for 0-angles away 
from 90° a grid of 100x100 starting points is sufficient 
to successfully simulate the resistivity. However, close to 
90° the orbits are rapidly changing with 9, and the in- 
terplane resistivity is dominated by a few small, closely- 
spaced orbits. In this case it is necessary to greatly in- 
crease the Fermi surface resolution, which in turn greatly 
lengthens the duration of the simulation. 

In performing these simulations the interest lies in their 
ability to reproduce the Lebed magic angle effect and the 
Yamaji oscillations, as it is these phenomena that need 
to be distinguished from one another. Less important 
are the Danner-Kang-Chaikin oscillations and the third 
angular effect. Although they too are reproduced by the 
Chambers formula, an analysis of these effects does not 
yield a great deal of useful information. That said, the 
in-plane peak effect, which is intimately related to both 
the Danner-Kang-Chaikin and third angular oscillations, 
is of great interest, but will be dealt with in a different 
manner in Section IVDI 

Simulations of the angle-dependence of the interplane 
resistivity at 42 T and several values of the azimuthal an- 
gle, (j>, for the Q1D sections of Fermi surface are shown in 
Figure El At — ±90°, which corresponds to the mag- 
netic field lying parallel to the Q1D sheets, the Lebed 
magic angle effect can be clearly seen as dips in the mag- 
netoresistance. As the </>-angle is changed the frequency 
of the dips also changes. At low (^-angles the amplitude 
of the dips drops, and at (j> — 0° they are no longer 
observed. The Danner-Chaikin oscillations are seen as 
smaller features near 9 = 90° at low azimuthal angles. 

The validity of these simulations can be checked by 
calculating the frequency in (tan^)" 1 , 1/x, of the Lebed 
magic angle dips for each value of the azimuthal angle. 
In this way the ^-dependence of \ can be fitted to the 
equation 



where (j>o corresponds to the magnetic field lying parallel 
to the Q1D sheets jU- The results of such an analysis are 
shown in Figure [7| From the fit the value of Xo: which is 
equal to the c lattice parameter divided by the interlayer 
distance, d±, is found to be 0.8658 ± 0.0005. This can be 
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FIG. 6: The simulated interplane resistance resulting from 
solving the Chambers formula numerically for the Q1D Fermi 
surface sheets of re-(ET)2Cu(NCS)2, as described by Equa- 
tion H3 The ^-dependences are shown for a fixed magnetic 
field of 42 T and a selection of values of the azimuthal angle, 
<f>. The Lebed magic angle effect dominates the magnetoresis- 
tance, except at low 0-angles where the Danner-Kang-Chaikin 
oscillations are seen around 8 = 90°. 

compared with the value of c/d± = 0.861±0.001 obtained 
from X-ray scattering measurements p|- 

The inset to Figure [7| shows the results of analysing 
the Lebed magic angles if they were mistakenly taken 
for Yamaji oscillations. The dotted line is a fit to the 
equation 

*r*W = [fc>s 2 (0-£) + ^in 2 (0-O] 1/2 , (10) 
which is valid for a Fermi surface with an elliptical cross- 
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FIG. 7: The (^-dependence of x, deduced from the frequency 
of the simulated Lebed magic angle dips. The dotted line is 
a fit to Equation |5] The insert shows the polar plot of fcjj" ax 
versus 4> that would result if the Q1D features were mistaken 
for Q2D Yamaji oscillations, with the dotted line representing 
a fit to Equation I10I 

section and where fc|j nax is the maximum in-plane Fermi 
wavevector projected on the plane of rotation of the field 
and is found from the frequency of the Yamaji oscilla- 
tions. k a and are the major and minor semi-axes of 
the Q2D Fermi surface pocket respectively j^. This fit 
suggests the existence of an elliptical Q2D pocket with 
major and minor axes of 2.37 nm _1 and 0.29 ran -1 re- 
spectively. If these were experimental results then it is 
easy to see that in the absence of any other evidence 
such a Fermi surface pocket might seem quite reasonable. 
However, the mistake becomes apparent when the funda- 
mental frequency of the quantum oscillations that would 
be expected from a closed pocket of this size, 227 T, is 
compared with the experimentally determined value of 
599 T. 

Figure [S] shows the simulated interplane magnetoresis- 
tance that arises from the Q2D closed Fermi surface pock- 
ets at 42 T and various ^-angles. It is seen that the traces 
are dominated by the peaks of the Yamaji oscillations. 
/c|| lax can be extracted from the frequency of the oscil- 
lations at each 0-angle, and the result of fitting this to 
Equation is shown in Figure The resulting Q2D 
pocket has a major axis, kf,, of 2.476 ± 0.001 nm _1 , and 
a minor axis, k c , of 0.733 ± 0.002 nm _1 . This would 
give rise to quantum oscillations with a fundamental fre- 
quency of 598 T, which is in agreement with the value 
measured from the SdH effect. 

It is now possible to mark the differences expected be- 
tween the shape of the Q2D pocket that results from the 
correct analysis of the Yamaji angles, and that from the 
mistaken identification of the Lebed magic angles. The 
obvious difference is that the pockets are perpendicular to 
each other, with the long axis of the true Q2D pocket ly- 
ing along the <fi = 0° direction. If the samples used in the 
experiments had been oriented by optical measurements 




(degrees) 

FIG. 8: The simulated interplane resistance resulting from 
solving the Chambers formula numerically for the Q2D Fermi 
surface pockets of re-(ET)2Cu(NCS)2, as described by Equa- 
tion H>] The ^-dependences are shown for a fixed magnetic 
field of 42 T and a selection of values of the azimuthal angle, 
<j>. The Yamaji oscillations dominate the magnetoresistance, 
taking their maximum frequency when the in-plane field is 
applied along <fi — 0°, which is parallel to the major axis of 
the Q2D pocket. 



then this would be sufficient to distinguish the AMROs. 
However, this is not the case. The major axes of the two 
alleged pockets are similar to each other, and an exper- 
imental error is likely to encompass them both. Thus it 
is to the minor axis that one must look to separate the 
two AMRO effects. 

All in all the simulations agree very well with the ex- 
perimental results of both X-ray scattering and the SdH 
effect. Nevertheless, it is worthwhile to look more closely 
at the various errors that might be introduced into the 
simulation process along the way. The first, most general 
problem to be addressed is that the simulations are semi- 
classical, and take no account of the quantum oscillations 
and, more importantly, the magnetic breakdown. As the 
perpendicular field is increased to high magnitudes, the 
experimentally measured AMRO will become affected by 
magnetic breakdown, as more and more Q1D carriers 
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FIG. 9: Left: The data points are the fc™ ax values as obtained 
from the frequency of the Q2D Yamaji oscillations in the sim- 
ulated resistance, the dotted line is a fit to Equation I1UI and 
the solid line is the resulting Fermi surface pocket. Right: A 
reminder of the in-plane Fermi surface of k-(ET)2Cu(NCS)2- 



tunnel through the energy gap and become Q2D carri- 
ers. Eventually the system will resemble one large Q2D 
Fermi surface pocket whose cross-section in the highly 
conducting planes is the /3-orbit. However, AMRO tend 
to be most concentrated near to 9 = 90°. In fact, the 
actual experiments were performed at 42 T, and it will 
be seen that almost all the important AMRO features 
occur at 6>-angles of around 70° or higher. The perpen- 
dicular magnetic field at B = 42 T, 9 = 70° is such that 
for the ti8-sample less than one quasiparticle in fourteen 
has sufficient energy to bridge the gap between the Fermi 
surfaces. For h8 this value is less than one in fifty, and 
the probability of breakdown for both types of sample 
decreases towards zero as 9 approaches 90° (whether or 
not the probability actually reaches zero at 0=90° de- 
pends on the relative sizes of t a and E g ). Thus for the 
current situation the magnetic breakdown is only a minor 
consideration. 

The most likely entry point for errors to make their 
way into the calculations is via the values chosen to rep- 
resent the various physical parameters. It has already 
been mentioned that the values chosen for tb, t c and Ep 
used in conjunction with the effective dimer model repro- 
duce the measured Fermi surface very well, so attention is 
turned to the other parameters, namely E g and t a . The 
value of 7.8 meV chosen for E g is derived from a mea- 
sured value of the magnetic breakdown, which has a large 
error associated with it 20]. However, a quick glance at 
how such an error propagates reveals it to be relatively 
unimportant: E g represents the gap in fc-space between 
the Q1D and Q2D sections of Fermi surface, and an order 
of magnitude estimate of this gap in terms of wavevector, 
Afc, is given by Afc/fcp ~ E g /Ep 0. Using estimates 
for kp and Ep 3], it is found that Afc ~ 0.4 nm -1 . As 
the area of the Q2D pocket is well defined, any error on 
the size of the gap would lead to errors in kb and fc c , the 
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FIG. 10: The angle-dependent interlayer magnetoresistance 
of h8 k-(ET)2Cu(NCS)2 at various values of the azimuthal 
angle, 6. T « 500 mK and B = 42 T. 



axes that define the pocket. A generous error on kb is 
±0.2 nm" 1 , or half of Afc, which represents what would 
happen if the energy gap were allowed to be zero. Fixing 
the area, this leads to a error of around ±0.1 nm" 1 on 
fc c . Even with such an uncertainty on the magnitude of 
the minor axis, it would still be possible to distinguish 
between the results arising from the Yamaji oscillations 
and those from the Lebed magic angles. 



The value used to represent t a , the transfer integral 
along the crystallographic a direction, is based on a pre- 
liminary analysis of the in-plane peak effect Q. As it is 
quite small, 0.04 meV, it is likely to have associated with 
it a significant relative error. However, the positions of 
the AMRO features arising from the Yamaji and Lebed 
effects are unaffected by the magnitude of the t a param- 
eter. In the current situation the amplitudes of the oscil- 
lations are of little concern, thus, for the moment, neither 
is the precise value t a . 
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FIG. 11: Left: The value of x> obtained from the frequency of the resistance dips caused by the Lebed magic angle effect, 
at various values of the calibrated azimuthal angle, (p. The dotted line is a fit to Equation Right: The values of fcj| lax (0) 
obtained from an Yamaji analysis of the peaks in resistance at various values of the calibrated azimuthal angle, <f). The dotted 
line is a fit to Equation 1101 constrained so that the resulting Q2D pocket (solid line) has an area corresponding to the measured 
fundamental frequency. In both figures the squares are the data from hS sample #1, the circles are the data from h8 sample 
#2 



2. Experiments 

Figure ^| shows a selection of the measured angle- 
dependences of h8 k-(ET) 2 Cu(NCS)2, at various values 
of the azimuthal angle, in a field of 42 T and at temper- 
atures around 500 mK. In order to analyse such angle- 
dependences, the position of each AMRO peak and dip 
is recorded. The frequency in (tan#) _1 of the peaks and 
dips at each jangle is then found for each sample, and 
the results are compared to those obtained from the sim- 
ulated resistance. The AMRO arising from the Q1D and 
Q2D Fermi surfaces are thus identified and the measured 
0-angle can be calibrated so that <f> — 0° is perpendicular 
to Q1D sheets. The frequencies of the dips arising from 
the Q1D Lebed magic angle effect for each sample are 
combined and fitted to Equation [5] The result is shown 
in the left hand side of Figure ^] The fit is good and it 
is found that xo{h8) = c/dj_ = 0.89 ± 0.10, which is in 
reasonable agreement with value of 0.861 ± 0.001 found 
from X-ray scattering Q. 

The results of calculating fc|j nax (0) for the resistance 
peaks arising from the Q2D Yamaji oscillations for each 
sample are also combined, and these data, together with 
the curve obtained by fitting to Equation 1101 are shown 
in the right hand side of Figure ^] 

It is seen that almost all the fc™ ax (0) data are con- 
centrated around the region where the magnetic field 
is roughly perpendicular to the flatfish portion of the 
Q2D pocket — at the other 0-angles the resistance is 
dominated by the Lebed effect. This means that the 
major axis, fcf,, is ill-defined and it is necessary when 
performing the fit to fix the area of the pocket so that 




FIG. 12: The result of analysing all the AMRO as if they 
were Yamaji oscillations. The squares are the data from h8 
sample #1, the circles are the data from h8 sample #2, the 
hollow symbols are data from Lebed magic angle dips, the 
solid symbols are data from Yamaji peaks, and the dotted 
lines are the fitted curves from Figure 1111 This shows that 
the Q1D effects dominate when the field is nearly perpendic- 
ular to the sheets, and the Q2D effects dominate when the 
field is nearly perpendicular to the flatter edge of the Fermi 
surface pockets. 



it reproduces the measured fundamental frequency of 
the SdH oscillations. It can be seen from the figure 
that this fit is reasonable; and the results obtained are 
k c (h8) = O.SOiO.OSnm- 1 and k b (h8) = 2.28±0.15 nm" 1 . 

It is illustrative to calculate /C|j nax (0) for all the h8 data 
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FIG. 13: ^-dependence of the simulated interlayer resistance 
at 6 = 90° for both the Q1D and Q2D sections of Fermi 
surface at a fixed field of 42 T. This is the same geometry as 
the third angular effect. 



(Lebed angles and Yamaji oscillations) and display the 
results on the same polar plot. In this way it is easy to 
see which Fermi surface section dominates the resistance 
at a given azimuthal angle. This plot is shown in Fig- 
ure El The hollow symbols are the Lebed magic angle 
dips, and the solid symbols are the Yamaji peaks. The 
dotted lines are the fitted curves from the previous two 
figures. It is clear that the Q1D Fermi surface sheets 
dominate the angle-dependent magnetoresistance when 
the field is roughly perpendicular to the sheets, and the 
Q2D pockets dominate when the field is roughly perpen- 
dicular to their flatter edges. 

In order to shed light on this behaviour, Figure ITTA 
shows the result of simulating the interlayer resistance at 
fixed field of 42 T and a 0-angle of 90° over the whole 
range of azimuthal angle, for both the Q1D and Q2D 
Fermi surface sections. This is just the simulated third 
angle effect plotted in polar coordinates, and is chosen 
to be representative of the magnitude of the resistance 
at 42 T in the angular region (70° -90°) over which 
AMRO are observed. By qualitatively comparing Fig- 
ures El an d El it is seen that the range of azimuthal 
angles over which the Q1D Fermi surface dominates the 
measured magnetoresistance is similar to that over which 
the simulated R zz (0 — 90°) due to the Q1D sheets is 
lower than that due to the Q2D pockets. The inverse is 
true for the range of angles where the measured magne- 
toresistance is dominated by Q2D AMRO effects. Given 
that the simulations have already been shown to be trust- 
worthy by reliably reproducing experimental data, it ap- 



TABLE I: Comparison of results derived from the AMRO 
of h8 and d8 k-(ET)2Cu(NCS)2, and from the semi-classical 
transport simulations. 





Xo 


k b (nm _i ) 


k c (nm L ) 


h8 


0.89 ±0.10 


2.28 ± 0.15 


0.80 ±0.05 


d8 


0.70 ±0.15 


2.19 ±0.35 


0.83 ±0.13 


simulations 


0.8658 ± 0.0005 


2.476 ± 0.001 


0.733 ± 0.002 



pears that the resistances from the two sections of Fermi 
surface combine in similar way to resistors in parallel, i.e. 
the overall resistance of the system at a given 0-angle is 
dominated by the section of Fermi surface that takes the 
lowest resistance at that angle. However, the situation 
away from 8 = 90° is not quite as simple as the parallel re- 
sistors scenario, as it is found that adding the simulated 
angle-dependent resistances for each Fermi surface sec- 
tion using i2j^ ta j = -Rqid + i?Q2D does not successfully 
reproduce the experimental results. 

Similar AMRO effects were measured in the d8 sam- 
ple. The left hand side of Figure ITU shows the azimuthal 
angle-dependence of the frequency of the Lebed magic 
angle dips. The fit to Equation is reasonable but not 
nearly as good as that for the h8 sample. The value of 
Xo{d8) = c/d± is found to be 0.70 ± 0.15. This can be 
compared to the value of 0.862 ± 0.001 found from X-ray 
scattering measurements of the deuterated salt j3j|. It 
is not entirely clear why the results are not as successful 
for the d8 sample as for the h8, however it is seen that 
there are nearly half the number of data points in the d8 
fit than the h8, and further, the density of points in the 
regions around <j) = 90°, where the xo parameter is best 
defined, is much lower in the case of the d8 fit. 

The right hand side of Figure El shows the azimuthal 
angle-dependence of k™ ax as calculated from the fre- 
quency of the Yamaji oscillations. In order to fit the 
data to Equation El the area is again constrained to 
produce quantum oscillations of the correct frequency. 
The axes of the pocket are thus found to be; k c (d8) = 
0.83 ± 0.13 nm^ 1 and k b (d8) = 2.19 ± 0.35 nm" 1 . 

A comparison of the results of analysing the AMRO 
measured in the h8 and d8 samples, and those simulated 
using the Chambers Formula, is shown in Table [I] It is 
seen that the results for h8 and d8 agree with each other 
to within the error ranges. It is also seen that there is 
reasonable correlation between the experimentally deter- 
mined values and those from the simulations. This im- 
plies that the parameters used in the simulation program 
are good approximations to the real values, and that it is 
possible to explain the AMRO in terms of purely semi- 
classical effects. 



D. Characterising the interlayer transport 

For a system with a three-dimensional Fermi surface, a 
series of quasiparticle orbits are possible in the presence 
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FIG. 14: Results for the d8 sample. Left: The value of x> obtained from the frequency of the resistance dips caused by 
the Lebed magic angle effect, at various values of the calibrated azimuthal angle, (j>. The dotted line is a fit to Equation 
Right: The values of fc™ ax (</!)) obtained from an Yamaji analysis of the peaks in resistance at various values of the calibrated 
azimuthal angle, (/>. The dotted line is a fit to Equation 1101 constrained so that the resulting Q2D pocket (solid line) has an 
area corresponding to the measured fundamental frequency. 





FIG. 15: Examples of the orbits possible on the Q2D (left) 
and Q1D (right) sections of the Fermi surface defined by 
Equation |H] when the magnetic field is applied parallel to 
the highly conducting planes. All the orbits tend to average 
the interlayer velocity to zero, and hence produce an increase 
in the interlayer resistance. For the purposes of the illus- 
tration the interlayer transfer integral has been exaggerated 
compared to its experimentally determined value. 



of an exactly in-plane magnetic field, many of which are 
very good at averaging the interlayer velocity towards 
zero (see Figure IT5|) . This is the origin of the in-plane 
peak effect in p zz , and s ugg ests a coherent nature to the 
interlayer transport [3> HH • 

For highly anisotropic materials, the angular width, 
2 A, of the in-plane peak, when measured in radi- 
ans, can be approximated by 2v™ ax /i>|| , where v± ax 
is the maximum of the out-of-plane component of the 
quasiparticle velocity and U|| is the in-plane compo- 
nent parallel to the plane of rotation of the magnetic 
field [9|. The ^-dependence of V\\ can be calculated for 



k-(ET) 2 Cu(NCS)2 using the dispersion relation in Equa- 
tion and the in-plane Fermi surface parameters dis- 
cussed in Section ITTT1 v± is given by h~ 1 dE/dk z , and 
so is a constant equal to 2i a acos(/3 — ir/2)/H. In 

this way 2A is calculated, with the value of t a left as 
the only adjustable parameter. This value is determined 
by comparing the results of the calculation with the ex- 
perimentally derived values for the width of the in-plane 
peak. This is the same method used in Reference 0, 
however in the present case the monoclinic structure of k- 
(ET)2Cu(NCS)2 is taken into account, and hence slightly 
different results are achieved. 

The result for the h8 sample is shown in the left 
hand side of Figure ^5] Here the points are the ex- 
perimental data for two h8 samples, and the solid lines 
are the results of the calculations obtained by setting 
t a (h8) — 0.065 ± 0.007 meV. The continuous curve arises 
from closed orbits on the Q2D FS pocket, which are pos- 
sible when the magnetic field is directed along any </>- 
angle. The closed loops correspond to orbits about the 
Q1D FS sheets, which are only possible when the field 
is directed along a limited range of <j). Away from this 
0-range the data follows the continuous curve fairly well 
and agrees with the ^-calibration found in the previous 
section. Around <j> = 0°, 180°, and 360° the width of the 
peak can be governed by any of the three sets of closed 
orbits possible; those on the broadly curved, convex re- 
gion of Q1D sheets; those on the pointed, concave region 
located at the Brillouin zone boundary; or those on the 
Q2D pocket. The way in which these orbits will combine 
to produce the in-plane peak is not entirely clear, but it 
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FIG. 16: Azimuthal angle dependence of the angular width, 2A, of the in-plane peak for h8 (left) and d8 (right) k- 
(ET)2Cu(NCS)2 at B — 42 T, T ~ 500 mK. The points are the experimentally measured widths and the solid lines are 
the simulated width with t a (h8) = 0.065 ± 0.007 meV and t a (d8) = 0.045 ± 0.005 meV. The continuous curve represents the 
contribution from the Q2D Fermi surface pockets, whereas the closed curves are those from the Q1D sheets. 



might be expected that the conductivities of each orbit 
sum to produce the total conductivity, as is the case in 
the Chambers formula. To the first approximation the 
interlayer resistivity is found by inverting the interlayer 
conductivity and so, 

1 111 

~ Ozz = Vzz 1 +Czz 2 +Vzz 3 + ■■■ ~ 1 1 1" ••• 

PzZ PzZ\ PZZ2 PzZ 3 

(11) 

i.e. the resistivity contributions from each orbit combine 
like resistors in parallel, and it is the path with the small- 
est resistivity that dominates the total resistance. Thus, 
in terms of the in-plane peak effect it is likely that the or- 
bits that are the least efficient at averaging the interlayer 
velocity towards zero will dominate the resistance. 

The results for the d8 sample are shown in the right 
hand side of Figure ^| The best agreement between cal- 
culations and experiment was found by setting t a (d8) = 
0.045 ± 0.005 meV. 



VI. CONCLUSIONS 

In summary, several physical properties have been 
measured for both hydrogenated and deuterated k- 
(ET) 2 Cu(NCS)2- No disparity has been found in the size 
and shape of the Fermi surfaces of the two isotopes, their 
effective masses, their scattering rates or their energy 



level structures. The only discernable difference found 
was in the interlayer transfer integral, which appeared 
lower for the deuterated salt. 

The size of the interlayer warping is determined by the 
transfer integral in this direction. An increased warping 
means that the Fermi surface will be less able to nest. 
Thus, if indeed the superconductivity in this material 
is aided by nestability @, IU, El then the higher transfer 
integral in the ft.8-salt would help to explain its lower 
superconducting transition temperature. 

Further, it has been shown that the measured angle- 
dependent magnetoresistance oscillations can be repro- 
duced via purely semi-classical, Boltzmann transport 
considerations. The observed peak in the resistance in 
the presence of a nearly in-plane magnetic field suggests 
that at low temperatures and ambient pressure the Fermi 
surface of k-(ET)2Cu(NCS)2 is a three-dimensional ob- 
ject that extends throughout reciprocal space. Thus it is 
not necessary to invoke non-Fermi liquid effects in order 
to describe the angle-dependent interlayer transport in 
this material. 
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